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The actinide cubic Laves compounds NpAk>, NpOs2, NpFe2, and PuFe2 have been examined 
by X-ray magnetic circular dichroism (XMCD) at the actinide Ma,6 absorption edges and Os 1/2,3 
absorption edges. They have the interesting feature that the An — An spacing is close to the so-called 
Hill limit so that substantial hybridization between the 5/ states on neighboring atoms is expected 
to occur. The XMCD experiments performed at the M.4,5 absorption edges of Np and Pu allow us 
to determine the spectroscopic branching ratio, which gives information on the coupling scheme in 
these materials. In all materials the intermediate coupling scheme is found appropriate. Comparison 
with the SQUID data for NpOs2 and PuFe2 allows a determination of the individual orbital and 
spin magnetic moments and the magnetic dipole contribution m m d- The resulting orbital and spin 
magnetic moments are in good agreement with earlier values determined by neutron diffraction, and 
the values of m m d are non-negligible, and close to those predicted for intermediate coupling. There is 
a comparatively large induced moment on the Os atom in NpOs2 such that the Os contribution to the 
total moment per formula unit is ~30% of the total. The spin and orbital moments at the Os site are 
parallel, in contrast to the anti-parallel configuration of Os impurities in 3d ferromagnetic transition 
metals. Calculations using the LDA+[/ technique are reported. The ab initio computed XMCD 
spectra show good agreement with experimental spectra for small values (0-leV) of the Hubbard U 
parameter, which underpins that 5/ electrons in these compounds are relatively delocalized. The 
calculations confirm the sign and magnitude of the experimentally determined induced magnetic 
moments on the Os site in NpOs2. A posteriori, by comparison of the theoretical and measured 
XMCD spectra, we can determine the most appropriate LSDA+ U variant and the value of U. 

PACS numbers: 75.25.-j,78.70.Dm, 71.20.Lp, 71.15.Mb 



I. INTRODUCTION 

The central question of the solid-state physics of the 
actinides is the extent and influence of the hybridization 
of the 5/ electrons. The vast panoply of properties, from 
localized to itinerant systems, from semiconductors to su- 
perconductors, is caused by variations of this property, as 
a consequence of the atomic configuration and environ- 
ment. The 5/ electrons have an extended wavefunction 
in real space so that they can often interact with 5/ states 
of neighboring atoms; more often, however, hybridization 
occurs between the 5/ states and other electron states of 
the same or neighboring atoms. For example, in metallic 
systems the 5/'s almost always interact with the conduc- 
tion states, which are made up of 6d and 7s states. In 
compounds with extended p or d states the 5/'s can have 
complicated interactions with such states, and a surpris- 
ing number of different phenomena are observed in the 
solid-state properties [lHj]. 

Many years ago Hill @] pointed out that there was a 
connection between the properties of actinide compounds 



depending on the inter-actinide spacing dAn in the mate- 
rial. For small dAn the materials tended to be super- 
conductors, whereas for large dAn the materials tend 
to order magnetically. These ideas addressed the hy- 
bridization between 5/ states on neighboring sites. In 
today's language we would recognize that for small dAn 
the bandwidth is large, favoring a suppression of mag- 
netism and possible superconductivity, whereas for large 
dAn the bandwidth is small and there is a tendency to 
localization of the 5/ states, leading to magnetic order. 
Despite the Hill criterion now being 40 years old, the 
ideas are still of interest. The critical spacing where these 
phenomena change was estimated by Hill as (^71^3. 3 A. 

It was recognized about this time that the actinide 
cubic Laves phases (C15 face-centered cubic structure 
of MgCu2 type) were well-packed structures with dAn 
close to the critical spacing. A systematic study of these 
materials 0-13 1 in polycrystalline form was performed at 
Argonnc National Laboratory in the 1970s. The mea- 
surements reported were electrical resistivity, magneti- 
zation, Mossbauer spectroscopy, and neutron scattering. 
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We shall draw on this body of work, but note that these 
techniques were unable to give any detailed information 
on the individual spin and orbital moments on the ac- 
tinide (or other) ions. Since it is important to know these 
values to compare with theory, which is now capable of 
treating the strong hybridization observed, we have revis- 
ited some of these materials with x-ray magnetic circular 
dichroism (XMCD). XMCD is element and shell specific 
and can obtain the individual spin and orbital moments 
in the structure on different atoms. 

In this paper we shall concentrate on the ferromagnetic 
systems, NpAl 2 , NpOs 2 , NpFc2, and PuFe2. Some of 
their properties already known are shown in Table HI We 
use for PuFe2 the results from a 1988 study of a single 
crystal with polarized neutrons 

Important conclusions from this body of work were: 

1. NpAl 2 is probably a localized system {d,An is larger, 
the moment is larger, and Tq higher, than in 
NpOs 2 ). 

2. NpOs 2 is probably an itinerant system, with a low 
moment and Tq and almost negligible magnetic 
anisotropy. 

3. In the case of both NpFe 2 and PuFe 2 the actinide 
5/ states are mostly itinerant. 

4. There is hybridization between the actinide 5/ and 
Fe 3d states as shown by the polarized-neutron 
study of PuFe 2 and partial quenching of the or- 
bital moment at the Pu site. 

Starting in 1989 Eriksson, Johansson, and Brooks pub- 
lished a sequence of theory papers 16l-ll9l| addressing par- 
ticularly the magnetic properties of these Laves phase ac- 
tinide materials. These band-structure calculations were 
performed with the local-density approximation (LDA) 
and the large spin-orbit interaction was included. The 
results gave the ordered moments, as well as their con- 
stituent orbital (^l) and spin (us) parts. Since the mag- 
netic form factor, as measured in neutron scattering, de- 
pends on ratio of the orbital to spin contributions (20| . 
they also predicted the magnetic form factors. In 1993 
the question of orbital polarization within the local-spin- 
density approximation (LSDA) was considered [2l| giv- 
ing, by this time, a fairly complete theory. 

Since that time, a new powerful magnetic character- 
ization tool has emerged; the X-ray magnetic circular 
dichroism (XMCD) technique. This synchrotron based 
technique allows one to quantitatively estimate the spin 
and orbital magnetic moments of the absorbing atoms. 
XMCD has been demonstrated to work successfully for 
the 3d, 4c? and 5d transition metals. Indeed, a set of 
sum-rules relate the experimental integrated XANES and 
XMCD spectra to the expectation value of ground state 
operators. More recently, this technique has been used 
extensively for a number of uranium compounds at the 



M 4 , 5 absorption edges (M 4 = 3.73keV; M 5 = 3.55kcV) 
where the dipole-allowed signal corresponds to the pro- 
motion of an electron from a core 3d shell to the partly 
filled 5/ band. Such experiments have been performed 
on ferro mag netic uranium compounds with the NaCl 
structure [22 -25 1, the Laves phase ferromagnet UFe 2 [26|, 



heavy-fermion materials where the moment is induced 
with a magnetic field 27. 28 1. UTA1 ferromagnets j29j . 
and even multilayers containing uranium [30L l3lj |. Up 
to now, only two XMCD studies have been reported on 
the M4.5 absorption edges of a transuranium material, 
NpNiGa 5 ,by Okane et al. [H and work at the ESRF on 
Np 2 Coi 7 0. 

There have also been a number of theory efforts calcu- 
lating the band structures and comparing the resulting 
XAS and XMCD signals with experiment [33-139]. Us- 
ing atomic theory, Thole and Carra4(j42| have derived 
sum-rules that relate integrated intensities of XAS and 
XMCD to the ground-state expectation values of stan- 
dard operators, such as the orbital (L z ) and spin (S z ) 
magnetic moments of the absorbing atom. Moreover, 
atomic multiplet theory has been used to calculate the 
atomic core-level spectra at the M45 absorption edges, 
which is associated with the transitions f n — >• d 9 f n+1 
((43|. [32l| ) . However, there was no attempt to calculate 
the spectral shape of both the XAS and XMCD for nep- 
tunium using band-structure calculations. 

XMCD can also be performed on the actinides at the 
iV.4,5 edges, where the transitions are from the Ad core 
shell to the partially filled 5/ shell and the transition en- 
ergies are less than 1 kcV; however, these experiments are 
performed in the soft x-ray range. Due to the large core- 
hole lifetime broadening at the N edges, it is more diffi- 
cult to record an XMCD signal than with the M edges. 
Furthermore, at these low excitation energy the experi- 
ments are extremely surface sensitive, so it is not certain 
that bulk-like properties will be measured. For a recent 
example see Okane et al. in Ref . [44| . 

Although the XMCD technique is powerful, there are 
a few caveats to recall. The first is that although the 
orbital moment can be estimated directly, the spin mo- 
ment cannot. Rather than (S z ), the quantity deduced is 
(S e * ' ) where this is the effective spin moment. From the 
sum rules [40l-42j. we find that for the M4 5 absorption 
edges, 



(S e ") = (S z ) + 3{T Z ) 



(1) 



where (T z ) is the magnetic dipole operator, which de- 
scribes correlations between the spin and position of each 
electron. (T z ) cannot be easily estimated experimentally, 
but it may be calculated assuming the electron count and 
the coupling of the electron states [jj, LS, or intermediate 
coupling (IC)] . Tabulated values are given in Ref. [43| • On 
the other hand, for actinide systems it is often these ex- 
act quantities (electron count and coupling scheme) that 



3 



TABLE I. Compounds discussed in this paper. a is the lattice parameter, d,An is the actinide-actinide interatomic distance 
(V3/4)a , Tc is the Curie temperature, fi a is the spontaneous ordered magnetic moment determined by magnetization studies, 
is the magnetic moment at the actinide site determined by neutron diffraction, is the magnetic moment determined 
at the Np site by assuming the hyperfine field is proportional to the moment [13, and fir is the magnetic moment at the T site, 
determined in these cases by neutron diffraction. Recall that the Hill criteria^ gives the critical value of c^in~3.3A. 



Compound 


a D (A) 




TciK) 


jU (/Us/f.U.) 




M^p(A*s) 




NpAl 2 


7.785 


3.37 


56(1) 


1.21(1) 


1.50(5) 


1.51(4) 




NpOs 2 


7.528 


3.26 


7.5(5) 


0.44(3) 


0.25(5) 


0.40(4) 




NpFe 2 


7.144 


3.09 


492(8) 


3.38(10) 


1.09(8) 


0.87(6) 


1.35(5) 


PuFe 2 


7.150 


3.10 


564(8) 


3.4(2) 


0.39(2) 




1.73(1) 



are unknown. We can determine (T z ), equivalently the 
magnetic dipole contribution m m d=-6(T z ), by knowing 
precisely the magnetization since the total moment fitot 
= fiL + /^Sj and combining this with the information 
from XMCD experiments. There is considerable interest 
in this quantity for the actinides as for localized uranium 
systems (T z ) appears to be large and close to the IC 
value [23J, whereas the value of (T z ) seems to be close 
to zero for itinerant systems possessing a small magnetic 
anisotropy, such as UFe2 [26| ■ 

Determining the value of (T z ) by comparing the signal 
obtained in XMCD and the magnetization obtained by 
SQUID magnetometry may appear simple, but with the 
photon energies at the actinide M (or especially with N) 
absorption edges the signal is sensitive to surface effects 
on the sample. Thus the SQUID is measuring a bulk 
signal, whereas the XMCD measures the signal from the 
first 200nm, or even less in the case of the N absorption 
edges, of the sample. As we discuss below, it appears 
that magnetic domains near to the surface are strongly 
constrained. 

The objective of the present work is to increase our 
knowledge of the AnX2 systems shown in Table I. In par- 
ticular, we wish to establish the orbital and spin moments 
experimentally, and determine the (T z ) for transuranium 
systems. Comparison with the values found from the lat- 
est electronic-structure theories tests our understanding 
of the physics of these compounds, the role of the 5/ elec- 
trons, as well as reveals which computational approaches 
provide the best description. 

The plan of the paper is in Sec. II to give the exper- 
imental details, in Sec. Ill we present the experimental 
results, in Sec. IV we report the theoretical calculations 
and their results, and a discussion is presented in Sec. V. 



II. EXPERIMENTAL DETAILS 

The actinide Laves phases used in this study are all 
stoichiometric compounds and relatively easy to prepare 
in polycrystalline form by standard techniques. 

The polycrystalline samples of NpAl2, NpFe2 and 



NpOs2 were prepared by arc melting stoichiometric 
amounts of high-purity elemental constituents Np (99.9% 
pure), Al (purity 5N), Fe(purity 5N) and Os (99.95% 
pure) on a water-cooled copper hearth under an Ar (pu- 
rity 6N) atmosphere. A Zr alloy was used as an oxy- 
gen getter. The samples were melted five times to im- 
prove the homogeneity and the mass loss was below 0.5%. 
Phase analysis of the obtained ingots was performed. 
Crystallographic analyses were performed at room tem- 
perature by x-ray diffraction (XRD) on samples with a 
mass of about 25mg, finely ground and dispersed on a Si 
wafer. Data were collected in back-reflection mode with a 
Bruker D8 diffractometer and Cu Ka radiation selected 
by a Ge(lll) monochromator. A one-dimensional po- 
sition sensitive detector was used to cover the angular 
range from 15 to 120° with incremental steps of 0.0085°. 
The analysis of the XRD pattern confirms a cubic C-15 
Laves phase structure. In all cases x-ray diffraction was 
used to identify that single-phase material was present, 
with lattice parameters consistent with those reported in 
Table I. In the case of PuFc2 , the sample had been pre- 
pared previously at ITU and pieces of that preparation 
were used for the SQUID and XMCD experiments. 

The NpAl.2 , NpFc2 and NpOs2 samples were then pol- 
ished to get platelets with two flat surfaces necessary 
for XMCD experiments and sealed in a specific capsule. 
The samples, which were between 10-20mg (except PuFc2 
which was around 0.2mg), were glued with epoxy resin 
in the middle of the Al support. A 100/im thick Be win- 
dow, covering the sample, was then glued with sty cast 
on the support. For safety reasons, the Be window was 
then covered with a 13/mi thin kapton foil. Finally, the 
capsule was closed with an Al upper part, screwed on the 
support. All the encapsulation procedure was done un- 
der helium atmosphere, to avoid any surface oxidation. 
Those capsules are sealed at ambient pressure and sus- 
tain liquid helium temperature. From the NpOs2 batch, 
we kept one sample piece as-grown, which we will label 
the as-cast NpOs2 sample. 

Magnetization and magnetic-susceptibility measure- 
ments were carried out in the temperature range between 
2 and 300 K and in magnetic fields up to 7 T using 
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a Quantum Design MPMS-7 superconducting-quantum- 
interfcrcnce-device (SQUID) magnetometer. 

The x-ray-absorption-spectroscopy (XAS) and XMCD 
experiments were carried out at the ID12 beamline of the 
European Synchrotron Radiation Facility (ESRF) , which 
is dedicated to polarization-dependent spectroscopy in 
the photon-energy range from 2 to 15 keV(45[. For the 
experiments at the M 4:5 absorption edges of Np and Pu 
(3.6-4.1 keV), the source was the helical undulator Helios- 
II, which provides a high flux of circularly-polarized x-ray 
photons with a polarization rate in excess of 0.95. After 
monochromatization with a double-crystal-Si(lll), the 
rate of circular polarization is reduced to about 0.42 at 
the M5 edge and 0.50 at the M4 absorption edges for Np 
and to about 0.46 at the M5 edge and 0.54 at the M4 
absorption edges for Pu. The x-ray-absorption spectra 
were recorded using the total-fluorescence-yield detection 
mode in backscattering geometry for parallel fj, + (E) and 
antiparallel [i~ (E) alignments of the photon helicity with 
respect to a 17T external magnetic field applied along the 
beam direction. Element selective magnetization curves 
were recorded by monitoring the intensity of the XMCD 
signal at a given photon energy, as a function of the ap- 
plied magnetic field. 

The x-ray-absorption spectra for right and left circu- 
larly polarized x-ray beams were done assuming semi- 
infinite samples, but taking into account the various 
background contributions (fluorescence of subshclls and 
matrix as well as coherent and incoherent scattering) , the 
angle of incidence of the x-ray beam, and finally the solid 
angle of the detector 46|- . 

The Np and Pu edge-jump intensity ratio M5/M4 (de- 
fined as the ratio between the occupation numbers for 
the two spin-orbit-split core levels j = 3/2 and 5/2) were 
then both normalized to 1.57 according to the theoret- 
ical edge-jump ratio tabulated in the XCOM tables by 
Berger et al. in Ref. 

The XMCD spectra ^ + {E) — n~{E) were obtained as 
the difference' of the corrected x-ray-absorption spectra. 
To make sure that the final XMCD spectra are free of ex- 
perimental artifacts, measurements were also performed 
for the opposite direction of the applied magnetic field. 
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FIG. 1. (Color online) Temperature dependence of the mag- 
netization measured on NpAl2, NpFe2 and NpOs2 polycrys- 
talline samples from 4 to 300 K in a field noH = 7 T. The 
inset shows the magnetic field dependence of the magnetiza- 
tion measured on NpAl2, NpFe2 and NpOs2 at 5 K. 



work as we shall compare magnetization and SQUID data 
at B = 7 T. It was already noted that NpOs2 has a large 
high-field susceptibility [111 ], For the PuFe2 sample, the 
ordering temperature and magnetization were the same 
as reported earlier 
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Absorption spectra and branching ratios at the 
actinide M edges 



Isotropic absorption spectra from the three NpX2 sam- 
ples are shown in Fig. [2] (top panel) with that for P11FC2 
in the lower panel. Since the samples all have cubic 
symmetry, we can compare the x-ray absorption signals 
recorded at the M^-edges directly. Clearly, the isotropic 
x-ray absorption spectra for all three NpX2 samples are 
similar. From the branching ratio, B = Im 5 / (Im 5 + 1 M 4 ), 
where Im 4 5 is the integrated intensity of the isotropic 
white lines at the M^s edge, we may determine the ex- 
pectation value of the angular part of the valence states 
spin-orbit operator, (i/j\£ ■ s\ip) = 3/2(W 110 ), as 5l| : 



III. EXPERIMENTAL RESULTS 
A. Macroscopic magnetic characterization 

Fig. [T] shows the temperature (T) and magnetic field 
(HqH) dependence of the magnetization obtained with 
the SQUID magnetometer for NpAl 2 , NpFe 2 , NpOs 2 . 
The ordering temperatures and magnetization are sim- 
ilar to those reported by Aldred et al. (jl-Hljl . Note in this 
respect that the value for Np compounds in Table U have 
been extrapolated to B = T, whereas no attempt has 
been made to perform that extrapolation in the present 



3n h 2 l 5^ 



A 



(2) 



where rih is the number of holes in the 5/ shell, and A 
is a quantity dependent on the electronic configuration. 
Free-ion values for A have been calculated [Hl|, and are 
-0.005 for 5/ 4 and zero for 5/ 5 . 



Mossbauer spectroscopy on the NpX 2 materials [llj has 
shown that the ionic state of Np is close to 5/ 4 , and the 
neutron experiments reported on PuFe2 [H[ have shown 
that the ionic state of Pu is close to 5/ 5 . The correspond- 
ing number of 5/ holes are therefore 10 for the NpX2 
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FIG. 2. (Color online) Isotropic x-ray absorption spectra 
recorded at the Ms and M4 Np edges for the polished NpX2 
samples (top panel (a)) and recorded at the Pu edges for the 
as-cast PuFe2 sample (lower panel (b)) as a function of inci- 
dent photon energy. The spectra have been corrected for self 
absorption effects. 



and 9 for PuFe2- Density functional theory calculations 
(see Sec. IV) give values of 10.3 and 9.1 for actinides 
in NpOs2 and PuFe2, respectively, close to the expected 
ionic states. Moreover, since (W 110 ) = n-j / 2 - 4/3 n 5 / 2 , 
where the number of electrons in the individual shells 
corresponding to j = I ± s, i.e. j = 7/2 and j = 5/2 are 
n 7 / 2 and n 5 / 2 , respectively, we may determine the oc- 
cupation of the individual spin-orbit split electron shells 
[H l5l| . These quantities are presented in Table HU 

From Tabic HH we can observe that the branching ra- 
tio B, provided by the XAS experiment, for Np in the 
various NpX2 compounds are all similar and within the 
error bar is equal to that for Np metal ones which is 0.74 
[3]. Further, no difference has been observed between 
the as-cast and polished NpOs2 samples. One obtains 
a 5/ spin-orbit interaction per hole of about -0.36, very 
close to the value calculated in the intermediate-coupling 



TABLE II. Experimental branching ratio B deduced from the 
x-ray absorption spectra, together with derived 5/ -electron 
contribution to the valence spin-orbit interaction per hole 
(W /110 )/n; 1 -A;the electronic 5/-state occupations n e and the 
electron occupation of the j — 5/2 and j = 7/2 sub-shells. 
Experimental error bars are ~ ±1% 



Compound 


B 


(W lw )/n h -A 


nV 


"5/2 


"7/2 


NpAl 2 


0.738 


-0.345 


4 


3.21 


0.79 


NpFe 2 


0.742 


-0.355 


4 


3.26 


0.74 


NpOs 2 


0.750 


-0.375 


4 


3.34 


0.66 


PuFe 2 


0.803 


-0.507 


5 


4.10 


0.90 



approximation for the 5/ 4 electronic configuration 51 
This supports the conclusion drawn from the Mossbauer 
isomer shift. In the case of PuFe2, the branching ratio is 
about 0.80 which is close to the Pu metal value of 0.82 
[l|. Furthermore, it agrees well with the value calculated 
in the intermediate-coupling approximation for the 5/ 5 
electronic configuration [5l|. We notice that going from 
Np to Pu the number of electrons in the 5/ 7 / 2 sub-shell 
remains approximately the same, implying that by going 
from Np to Pu, we are filling with 1 electron more the 
5/5/2 sub-shell. 



C. Element specific magnetization curve recorded by 
XMCD 



In Fig. [3] are shown the XMCD signals recorded at 
the Np M4 : 5 absorption edges under magnetic field of 6T 
and at 10K in NpX2 compounds. We observe that the 
XMCD signal at the Np M4 edge is large and that its 
spectral shape consists of a single nearly symmetric neg- 
ative peak that has no distinct structure. This dichroic 
signal at the M4 edge is also characteristic for all ura- 
nium systems. The dichroic signal at the M5 edge is 
nearly three times smaller than at the M4. The spectral 
shape of the XMCD signal has an asymmetric s shape 
with two peaks - a negative and a positive peak. How- 
ever this positive peak is most pronounced in NpFe2 and 
quasi absent for NpOs2, as found in NpNiGas [32J. This 
asymmetric spectral shape, like for the M 5 -edge of U, de- 
pends strongly on the hybridization, Coulomb, exchange, 
and crystal-field interactions. We will see in section IV 
that within a band structure model, this asymmetry can 
be well explained [39| . 

Since one of the principal aims of the present investi- 
gation is to compare the moment determined by SQUID 
magnetization (Fig. [I} with that deduced by the XMCD 
technique, we show in Fig. [4] the XMCD signal at the M 4 
absorption edge as a function of applied magnetic field 
from polished samples of the NpX2 systems. 

It is clear in comparing Fig. 0]with Fig. [T]that the el- 
ement specific magnetization curve recorded at the max- 
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FIG. 3. (Color online) XMCD spectra measured at the Np M 5 
and M4 absorption edges for the polished NpX2 samples in an 
applied field of 6T at 10K. The spectra have been corrected 
for self absorption effects and for the incomplete circular po- 
larization rate. 



imum XMCD signal at the M4 absorption edge has a 
quite different form, and shows no sign of saturation for 
all the three NpAl2, NpFe2 and NpOs2 polished samples. 
Surprisingly, one observes opening of the hysteresis loops. 
It seems, therefore, that magnetic domains at the sample 
surface (probing depth at the M4 edge resonance is in the 
range of ^200nm) are strongly constrained, and do not 
rotate with the applied field. We have then recorded the 
element-specific magnetization curve at the maximum 
XMCD signal recorded at Fe K edge and Os L3 edge. 
Since these signals are taken with photons of higher en- 
ergy, the probing depth is larger, especially at the Os L 3 
edge. They are clearly closer to saturation, as would be 
expected, as they penetrate further (at least one micron). 
It also interesting to notice that despite the unexpected 
shape of the hysteresis curves in Fig. 21 they do con- 
firm one of the points made in the earlier studies. The 
magnetic anisotropy of NpFe2 is very large, that of NpAl2 
smaller, and that of NpOs2 almost negligible. This is con- 
sistent with the strong coercivity shown by the expanded 
hysteresis loop in NpFc2 and its absence in the case of 
NpOs2. Further, the coercive field of the hysteresis curve 
recorded at the Fe K-edge, as well as at the Os L3-edge, 
is smaller than the one recorded at the Np M4-edge. It 
means that the strong anisotropy of the magnetization 
at the surface, most probably caused by the polishing, 
clearly penetrates into the sample, despite the coercive 
field being smaller. As a direct consequence, we cannot 
extract meaningful parameters (such as the spin, orbital 
magnetic moments, and (T z )) from these data, and make 
a comparison with macroscopic measurements. Because 
of the very high anisotropics known to exist for NpFc2 
and, to a lesser extent, NpAl2, (see Fig. 0]), we judged 



-1.5 
0.6 

0.4 

0.2 

0.0 

-0.2 

-0.4 

-0.6 
0.6 



(a) 

• 


' r 
• 

• 


• 

™," 


• 

NpAI 2 


,-•/ 
*• f 




/> / — n p ; 
7 / ■ 


/ // 

— 


i •' 
• 

• 

NpFe 2 


- 1 1 1 1 1 1 1 1 1 1 1 


i i i i i i i i i i i — 


=f 



10K 


NpOs 2 



0.006 
0.004 
0.002 X 

s 

o 

0.000 ° 

-n 

o 

-0.002 J 
a 
ca 

-0.004 " 

-0.006 
0.03 



-6 -5 -4 -3 -2 -1 1 2 3 4 5 6 

Magnetic Field (T) 

FIG. 4. (Color online) Element specific XMCD magneti- 
zation recorded at the maximum XMCD signal at the M4 
absorption edge (E=3846eV) measured at 10K for the pol- 
ished (a) NpAl 2 ,(b) NpFe 2 and (c) NpOs 2 samples. To- 
gether is also shown the element specific XMCD magneti- 
zation recorded at the maximum XMCD signal at the Fe K 
absorption edge (E=7I13.7eV) and at the Os L3 absorption 
edge (E=10877eV) for the polished NpFe 2 and NpOs 2 sam- 
ples respectively. 



it difficult to remove all the surface effects in these two 
materials, and concentrated on NpOs2. 

Let us now return to the element-specific magnetic 
characterization by XMCD. As it can be seen from Fig. 
[5j the amplitude of the XMCD signal as a function of the 
magnetic field recorded at the M4-edge is weaker for the 
polished sample of NpOs2 than for the as-cast sample. 
Thus, the polished surface affects not only the approach 
to saturation, but also the absolute value of the signal. 
If we superimpose now these element-specific magneti- 
zation curves for the as-cast sample with the SQUID 
measurements, excellent agreement is obtained over the 
whole range of applied field. 
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FIG. 5. (Color online) Element specific XMCD magnetiza- 
tion recorded at the maximum XMCD signal at the Np M 4 
absorption edge (E=3846eV) measured at 5K for the as-cast 
and polished NpOs2 samples. The earlier results from Aldred 
et al. [ll|] are from magnetization experiments. 



Regarding the sample of PuFe2, which was as-cast, (it 
is almost a single crystal since the sample came from the 
same boule as used for the single crystal studied by Wulff 
et al. in Ref. [3|)i the element specific magnetization 
curve recorded at the Pu M4-edge shows saturation and 
similar hysteresis as the macroscopic measurements. 

The conclusion from this XANES and XMCD charac- 
terization is that, for actinides systems, great care has to 
be taken to assure that surface effects do not play a role in 
comparing XMCD and magnetization results, especially 
of highly anisotropic materials. 



D. X-ray magnetic circular dichroism at An M 4 ,5 
absorption edges 

The orbital contribution to the 5/ magnetic moment 
[40| is obtained from the dichroic signal integrated over 
the pair of spin-orbit split excitations, A/m 5 + A7m 4 nor- 
malized to the isotropic x-ray-absorption spectrum, 



2n h 



l Ah 



-{AI Ah +AI Mi ) 



(3) 



A second sum rule (41 , 42 1 correlates a linear combina- 
tion of the partial dichroic signals AIm 5 and AJm 4 with 
the effective spin polarization (S e ^), which, in turn, is 
related to the spin operator through Eq. [TJ and given by 



(S z )+3(T z ) = j 



M 5 



I Mi 



-(A/m 6 



:AJ. 



Mi 



(4) 



where T z is the component along the quantization axis 
of the magnctic-dipole operator. This operator, corre- 
lating the spin and position of individual electrons, is 



associated with the asphericity of the electronic cloud, 
distorted by crystal-field or spin-orbit effects, and with 
the spin anisotropy. (T z ) is therefore correlated to the 



charge and magnetic anisotropy [22 1 . 

The orbital and spin components of the total magnetic 
moment [i = —((L z ) + 2{S Z )) [1b can then be obtained 
from XMCD spectra, together with an estimate of (T z ), 
if the value of the total moment [i and the occupation 
number of the 5/ shell are known. Note that the magne- 
tization data will give the total moment per formula unit, 
fJ-tot = MAn + 2yUx, where [ix is the contribution from the 
moments at the X sites, and fj,A n = f^/n + t^cond where 
fJ'cond, is the contribution, usually small, from the 6d7s 
conduction band. There have been several attempts to 
study ficond combining magnetization and neutron data 
[52| and magnetic Compton scattering with XMCD [25| , 
and a rough estimate is that it is about —10% of the 
total moment. The reason for the negative sign is that 
the conduction electrons are polarized by the spin con- 
tribution at the actinide site and, as is well known, this 
spin contribution is antiparallel to the total moment in 
the light actinides, as the moment is dominated by the 
large orbital moment 0] . Thus for NpOs2 , in making the 
comparison of the XMCD results with those from the 
magnetization, we shall assume that [i CO nd is -0.1 fipfp, 
so that /j, 5 J p = 1.1/ijVp will be compared to the signal de- 
duced from XMCD. In the case of PuFe2 we use the 5/ 
moment deduced from the polarized-neutron study of a 
single crystal [14 1. 

As discussed earlier, the total signal as deduced from 
XMCD is unreliable for NpAl2 and NpFe2 since there is 
no saturation of the signal (see Fig. [4]), so the deter- 
mination of spin and orbital moments as well as the T z 
contribution is restricted to NpOs2 and PuFe2 samples. 

The XANES and XMCD signal recorded at the M 4i5 
edges for the 2 compounds discussed in this section are 
shown in Figs. [HI and [7] 

We have first determined the 5/ orbital magnetic con- 
tribution for the NpOs2 and PuFe2 samples. The results 
are given in the Table IIIII As for uranium compounds, 
the 5/ orbital contribution is large and parallel to the ex- 
ternal magnetic field. The difficulty in determining the 
5/ spin magnetic moments, as already mentioned above, 
lies in the fact that the magnetic dipole contribution can- 
not be measured directly. The spin magnetic moment 
can be individually determined cither by taking the (T z ) 
contribution from DFT or multiplet calculation (for a 
free-ion without crystal-field interaction) or by combin- 
ing results from neutron diffraction and/or macroscopic 
measurements. Concentrating on NpOs2, for which we 
have XMCD data at both the Np M and Os L edges, 
the XMCD results, assuming rih = 10, give a 5 /-orbital 



magnetic moment /i L J 



5f-Np _ 



IM/ib/Np from Eq. (3). 



With a total magnetic moment fi 5 ^ Np = 0.46/is (tak- 
ing into account the Os magnetic contribution extracted 
from XMCD, see section III.E, and the spd estimated 
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TABLE III. Experimental results for the Np and Pu 5/orbital 
and spin magnetic moments deduced from XMCD for Np0s2 
(scaled to 7T) and PuFe2 as well as the m m d=-6{T z )fiB con- 
tribution deduced combining results extracted from XMCD 
and SQUID or polarized neutrons measurements compared 
together with theoretical expectation values. Note that here 
we compare the measured values to model atomic calcula- 
tions. These values will be compared to DFT claculations in 
Sec. IV. 





NpOs 2 


PuFe 2 


An (hb /An atom) 


1.09(5) 


1.98(5) 


M s {//" (Ms /An atom) 


-1.15(5) 


-1.96(5) 


5f— An / 5f — An 

Ml /M S e// 


—0.95(8) 


— 1.01(8) 


Mtot (^s/f-U.) 


0.60(1) 




nl s o r An (Ms /An atom) 


0.46(5) 
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Ms /_An (/W An atom) 


-0.63(7) 
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Ml /Ms 


-1.73(10) 


-1.25(10) 


Theory IC 


-1.90 


-1.35 


m md (fie/An atom) 
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FIG. 6. (Color online) XANES measured with right and left 
circularly polarized X-rays and associated XMCD spectrum 
recorded at the Np Ma,z absorption edges under 17T and at 
5K for the as-cast NpOs2 sample. The spectra have been cor- 
rected for self-absorption effects and for incomplete circular 
polarization rate. 



contributions of Np), we obtain a spin magnetic mo- 
ment on Np sites ^g~ Np = — 0.63/is/A^p. It is note- 
worthy that [i 5 /~ Np '/ 'n°/~ Np = —1.73 is close to the 
Np 3+ free-ion value (—1.90). From Eq. 01 we can then 
extract (T z ), which is +0.09. In fact, it is better to 
consider the magnetic dipole contribution m m d normal- 
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FIG. 7. (Color online) XANES measured with right and left 
circularly polarized X-rays and associated XMCD spectrum 
recorded at the Pu M± >5 absorption edges under 17T and at 
5K for the as-cast PuFe2 sample. The spectra have been cor- 
rected for self-absorption effects and for incomplete circular 
polarization rate. In the insert is shown the element specific 
XMCD magnetization recorded at the maximum XMCD sig- 
nal at the Pu M4, absorption edge (E=3846eV) measured at 
5K. 



ized by the spin moment (j^f on the Np atom, and this 
can be compared directly to calculated values in Table 
IIII1 m m d/fJ,g =+0.83 whereas in NP2C017 this value is 
+1.36. The closest theoretical value is that of the in- 
termediate coupling of +0.56. A similar analysis can 
be done for PuFe2, but instead of using macroscopic 
measurements, since we do not have access to the Fe 



magnetic moments, we used the results of Ref. [14l ] us 
ing neutron diffraction. From the XMCD results, as 
suming rih = 9, one obtains a 5/-orbital magnetic mo- 
ment L4^ J 



5f-Pu = 19SflB / Pu from Eq [3J The deduced 



spin moment on Pu sites fi s Pu = — 1.59/is/Pu, and 



the ratio fi£ 

,5/ 



5f-Np l£f-Np 



-1.25. Then for the Pu 



m m dj fJ- s J =+0.37, which is small and closest again to the 
IC value of -0.22. 

The neutron result for PuFc 2 for fi 5 /^ Pu /fi 5 s f ~ Pu is 
-1.20 ± 0.05, in excellent agreement with the XMCD 
result of —1.25 ± 0.05. These values are slightly be- 
low the IC theory value of -1.35 indicating a hybridiza- 
tion with the Fe 3d electrons, as already discussed in 
Ref. 14|. Moreover, the individual values for the orbital 
and spin moments determined by neutron diffraction are 
^5/-Pn = 2.3±o.3l4 B and f/ s f ~ Pu = -1.9 + 0. 3fi B , which 
are in reasonable agreement with the values determined 
from the XMCD in the table. We expect the neutron re- 
sults to be reliable for the ratio fi°/~ An j l^ 5 /~ An but not 
so accurate for the individual values. The XMCD tech- 
nique has better precision, but shows that the neutron 
values are reliable. 

The 5/ moment at the Np site in NpOs2 is 0.46 ± 
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0.03/is in good agreement with the value of 0.40 ± 
0.04/is deduced from the hypcrfine field measured in the 
Mossbauer experiments of Aldred et al. 

From the XMCD measurements, we can conclude that, 
similarly to uranium compounds, the largest contribu- 
tion to the 5/ magnetic moment arises from the orbital 
contribution, which is parallel to the external field. The 
spin contribution is antiparallcl to the orbital one. In 
both NpOs2 and PuFc2, the expectation value of the 
magnetic-dipolc contribution is small and is compatible 
with the value predicted for 5/ electrons calculated in 
intermediate-coupling (IC) approximations. 

NpOs2 has one of the smallest ordered moments of any 
Np-intcrmctallic compound, and this value has been dif- 
ficult for theory to reproduce. 



E. X-ray magnetic circular dichroism at Os Z/2,3 
absorption edges 



Figure M shows the XANES and XMCD signals 
recorded at the £2,3 absorption edges of Os for the as- 
cast NpOs2 polycrystalline sample, after corrections of 
self-absorption effects, the latter being small. 

Using a similar analysis procedure as described in Rcf. 
[53| . the resulting application of the sum rules at the 
L2.3 absorption edges of Os, using a number of 5d holes 
of 3.07[54j and neglecting the (T z ) contribution, gives for 
the Os 5d states a spin and orbital induced magnetic 
moment contributions under a field of 17T and 5K. We 
then re-scale the values to 7T applied field using the ele- 
ment selective magnetization curve shown in FigjBjD and 
obtain: 



= +0.013 ±0.005/x B /Os 



Ml 
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2.0 



= 1.5 
3 



<5 LO- 



OT 
UJ 



< 

X 



0.5 



0.0 



7/r 



Os L. 



Os L. 




(a) 



1 — 1 — 1 — ' — ih — > — 1 — 1 — 1 — > — r 

10840 10880 10920 12360 12400 12440 12480 



Fe„Co V :Os 

49 49 2 



0.3 
0.2 



-0.1 



0.0 
I- -0.1 



O 
D 



Photon energy (eV) 



0.04 
0.03 
0.02 
0.01 
0.00 
-0.01 
-0.02 
-0.03 
-0.04 



i 5K 




1 1 1 1 1 1 


(b) ■ 

1 1 1 1 1 1 



1.2 
0.8 
0.4 
0.0 
-0.4 
-0.8 
-1.2 



O 
O 



o 

Q. 

ID 
fl> 



-15 -10 -5 5 

Magnetic Field (T) 



10 



15 



FIG. 8. (Color online) (a) Isotropic XANES and XMCD spec- 
tra recorded at the Os Z/2,3 absorption edges under 17T and at 
5K for the as-cast NpOs2 sample and under 6T and 300K for 
the Co 4 gFe49V2 sample doped with 3% Os. The XMCD spec- 
tra recorded on the as-cast NpOs2 have been multiplied by a 
factor 3. In (b) is shown the element specific XMCD mag- 
netization recorded at the maximum XMCD signal at the Os 
L3 absorption edge (E=10877eV) (solid points) measured at 
5K together with the element specific XMCD magnetization 
recorded at the maximum XMCD signal at the Np M4 absorp- 
tion edge (E=3846eV) (open points) for the as-cast NpOs2 
sample. 



IV. ELECTRONIC STRUCTURE 
CALCULATIONS 



The Os 5<i contribution to the total magnetization is 
0.18/0.60 ~ 30%. The total moment (neglecting the sp 
contributions) is parallel to the Np moment. The in- 
duced spin contribution dominates over the orbital con- 
tribution, however, it is noteworthy that the induced or- 
bital and spin magnetic moments are aligned parallel. 
This is different from Os impurities in Permendur alloy 
(Fe4gCo4gV2), where the induced spin and orbital mag- 
netic moments are aligned antiparallel. Theoretically, 
Tyer et al, J5f|, show that for 5% Os (which has ~4 
5d holes in the atomic state) diluted in a Fe matrix, the 
spin moment should be positive, but the orbital moment 
should be negative, as we have found experimentally 



We have performed electronic structure calculations of 
the NpX2 (X=A1, Fe, Os) and PuFe2 compounds using 
the WIEN2k [56| package, which implements the full- 
potential linearized augmented plane waves method with 
local orbitals to solve the Kohn-Sham equations. To deal 
with the strongly correlated nature of unfilled /-electron 
shells, we have performed both local spin density approx- 
imation (LSDA) calculations as well as LSDA+£/ calcula- 
tions, which include explicitly the strong onsite Coulomb 
interactions in the 5 /-electron subsystem. We have 
tested both the most common variants of the LSDA+J7, 
namely the around mean-field (AMF [H?}) and fully local- 
ized limit (FLL [58|) formulations. The spin orbital (SO) 
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FIG. 9. (Color online) Spin- and atom-resolved density 
of states calculated for NpOs2 with the AMF variant of 
LSDA+t/ and U = l.OeV. 
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FIG. 10. (Color online) Spin- and atom-resolved density of 
states calculated for PuFe2 with the FLL variant of LSDA+!7 
and U = l.OeV. 



interaction, which is essential for actinide atoms, was in- 
cluded in the second variational step [59[ , where also the 
orbital potential was introduced, including its spin cross- 
term. In some of the calculations performed for NpOs2, 
see below, the SO interaction was explicitly switched off 
on selected atoms, to probe its influence on the spin and 
orbital moments of the two elements. The structure pa- 
rameters were set to 7.189A, 7.528A, 7.144A, and 7.785A 
for PuFe 2 , NpOs 2 , NpFe 2 , and NpAl 2 , respectively. The 
atomic sphere radii were set to 2.9 a.u., 2.5 a.u., 2.35 
a.u., and 2.5 a.u., for actinide and Os, Fe, and Al, re- 
spectively. The RK max parameter was set to 8.0, which 
leads to a basis size of approximately 500-550, depend- 
ing on lattice parameter and muffin-tin radii. In total, 
10000 k-points were used, leading to 726 k-points in the 
irreducible wedge of the Brillouin zone. These param- 
eters were checked to provide well-converged electronic 
structure calculations. 

In Figs. [S] and [TU1 we show the density of states 



(DOS) of NpOs 2 and PuFe 2 , respectively, for values of 
U and variants LSDA+C/; which provide the best agree- 
ment with experimental XMCD spectra, see below. One 
can observe large exchange splitting in the Np 5/ states, 
which is mainly responsible for the different spin-up and 
spin-down projection of the total DOS. The two spin pro- 
jected DOSes for Os atoms look very similar to each other 
reflecting the small induced spin moment on this atom 
in NpOs 2 , see Table llVl The unoccupied states are dom- 
inated by sharp peaks of Np-5/ character. It is transi- 
tions to these states that dominates the x-ray absorption 
spectra. A similar situation for unoccupied states applies 
also for the case of PuFc 2 , sec Fig. [10l The strongest fea- 
ture in the unoccupied states is formed by a 2eV wide 
Pu-5/ band starting approximately 2eV above the Fermi 
level. For the occupied states it is the exchange-split Fe- 
3d band that represents a dominant contribution. 

The XAS have been simulated by the initial state ap- 
proximation in which the absorption strength at the 
edges is a function of the projected density of states of un- 
occupied 5/ states of the actinide element. The energy- 
dependent dipolc transition matrix elements have been 
included. For a description of the implementation in the 



Figure [TT] sum- 



FLAPW basis we refer to Rcfs. 
marizcs the resulting XAS spectra of NpOs 2 and their 
difference - the XMCD. For every system treated we have 
performed several independent simulations - varying the 
strength of the onsite Coulomb repulsion characterized 
by the U parameter in the LSDA+C/ method. Typically 
we have calculated the spectra for U within the range 
from to 4 eV. The value of J was set to 0.6 eV for both 
Np and Pu 5/ electrons. 

Due to of the strong influence of the spin-orbit inter- 
action on the 3d core shell, the 3d 3 / 2 and 3d 5 / 2 core lev- 
els are split by approximately 200 eV. To calculate both 
edges within the same energy window would thus require 
evaluation of the unoccupied states to more than 200 eV 
beyond the Fermi level, which is not feasible in electronic 
structure codes using linearized basis sets. On the other 
hand, it is reasonable to assume that the unoccupied 5/ 
states will be almost entirely present within a few eV 
above the Fermi level (as is supported by DOS curves, 
Figs. [9] and [T0|) and that transitions to dipole- allowed 
continuum states would form a very wide feature-less 
contribution that would not have any visible influence 
on the spectra. A similar argument applies for the unoc- 
cupied p-states well above the Fermi level. Therefore, we 
have calculated each edge separately in its energy range, 
considering states up to approximately 50 eV above the 
Fermi level. A sudden monotonous decrease in the calcu- 
lated intensity approximately 50 eV above the edge onset 
is artificial and originates from this cut-off. Also a cer- 
tain inaccuracy can be expected at energies well beyond 
approximately 20 eV above the Fermi level due to the 
linearization errors present in the basis functions, some- 
thing which should be kept in mind, when interpreting 
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FIG. 11. (Color online) Theoretical and experimental XAS 
and XMCD of iV/4,5 absorption edges of NpOs2. Theoretical 
results are presented for various values of the U parameter in 
the LSDA+C/ potential, as indicated in the legend. 



the calculated spectral features far above the edge onset. 

All spectra were broadened with a Lorentzian of 2.5cV 
half-width at half-maximum, which provided good match 
with experimental widths of absorption peaks. This 
broadening can be interpreted as a combined effect of in- 
strumental broadening and life-time broadening, respec- 
tively. The intensity of the M 5 XAS edge was normalized 
to the same peak value as the experimental one. The 
same scaling factors were applied also for M\ edges and 
XMCD spectra, in order to preserve their relative mag- 
nitudes. 

Comparing the theoretical and experimental XAS 
spectra of NpOs2, the first-principles calculations repro- 
duce the experimental findings with good qualitative ac- 
curacy with regard to the peak shapes, distance and 
branching ratio. These features are only very weakly 
sensitive to the chosen model of treatment of the 5/ elec- 
trons. A certain degree of sensitivity on the strength of 
the onsite Coulomb correlations £7 can be seen on the in- 
tensity of the Mi peak. Best agreement between theory 
and experiment is obtained for small values of £7. 

The XMCD spectra reveal much more details. Gener- 
ally, the XMCD signal on the M 5 edge is very weak, com- 
pared to the XAS intensity. In the experimental spec- 
tra there has an s-shape, starting with negative values, 
switching the sign to a positive peak and monotonously 
decreasing to zero. The origin of such an XMCD shape 
has been outlined previously |39| . 

Considering the sum rules expression (Eqs. [3] and [4}, 
the sine-like shape implies that the energy-integral of the 
XMCD spectrum over the Ms edge will yield a very low 
value and, hence, it is mainly the M4 edge signal that will 
determine the local spin and orbital moment magnitudes. 
Indeed, the XMCD signal at the M4 edge is substantially 



larger, showing only a simple negative peak of rather 
large relative magnitude compared to the XAS M4 peak. 
This indicates that for the three Np-based compounds, 
the ratio of the orbital and spin magnetic moments will 
be approximately the same if (T z ) is not changed. This 
finding should remain valid also for NpAl 2 and NpFe2 de- 
spite that in experiments we have not reached magnetic 
saturation. 

The theoretical XMCD spectra at the M5 edge show an 
interesting, rather strong sensitivity on the chosen model 
of the description of the 5/ states. The small s-like sig- 
nal is best reproduced by low £7 values. Even the LSDA 
reproduces the sine-like XMCD signal, however, it over- 
estimates its magnitude. The sensitivity of the XMCD 
spectrum shape on the chosen 5/ model is much reduced 
at the M4 edge, although the magnitude changes visibly, 
when adding the onsite Coulomb correlation effects into 
the calculation. Our first-principles calculations thus re- 
veal a strong sensitivity of the computed XMCD spectra 
to the degree of 5/ electron localization, as described 
by the Coulomb £7. In comparison with an experimen- 
tal XMCD spectrum, this feature provides a means of 
establishing the appropriate Hubbard £7 value for a spe- 
cific actinide compound. Here we find that, generally, for 
the NpX 2 compounds, they appear to be best described 
by low £7 values (£7=0-1 eV), i.e. the 5f electrons behave 
mostly as itinerant. 

For the NpOs2 compound we have performed a more 
detailed study of the electronic structure, in order to 
compare to experimental findings, including the Os local 
magnetic properties, as extracted from the experiment. 
Namely, we have focused on the influence of the spin- 
orbital interaction of both Np and Os atoms on their spin 
and orbital magnetic moments. The WIEN2k package 
offers the possibility to selectively turn the spin-orbital 
interaction on and off within the atomic spheres of indi- 
vidual atoms. Using this feature, we have calculated the 
electronic structure properties of NpOs2 in four configu- 
rations: 1) full spin-orbital calculation, 2) no spin-orbital 
interaction, and 3) and 4) spin-orbital interaction present 
only on the Np or on the Os atom, respectively. 

Generally, using the AMF variant of the LSDA+£7 
causes the spin moment to decrease and the orbital mo- 
ment behaves non-monotonically. It grows quickly al- 
ready for small £7 and then decreases in magnitude. In 
the FLL variant of the LSDA+£7, both spin and orbital 
moments are less sensitive to £7, although there is a weak 
increase of the orbital moment magnitudes with growing 
£7. The (T z ) term expectation value is substantially in- 
fluenced by inclusion of the onsite Coulomb correlations, 
but the magnitude remains within the 0.1-0.2 range (on 
the Np atom). 

The Os spin magnetization follows the evolution of the 
spin magnetic moment on the Np atom, which indicates 
that the Os magnetization is of induced nature, far from 
the atomic regime and Hund's rules. The orbital mo- 
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TABLE IV. Summary of the magnetic characteristics of NpOs2 based on electronic structure calculations. The spin, orbital 
and total magnetic moments ([is, Mio Mtot) and magnetic dipole term (T z ), all in Bohr magnetons, are shown as a function of 
Coulomb U parameter and variant of the LSDA+(7 orbital functional, for both Np and Os atoms, respectively. The exchange 
J parameter is fixed at 0.6eV in all calculations (of course, except for the LSDA ones). We chose a sign convention in which the 
total moments per formula unit (pitot~ Np + 2/^| d ~° s ) are always positive. Recall from Table ITTT1 that the experimental values for 
Np 5/ states at 7T are hl = +1.09(5); us = -0.63(7); giving fi to t = +0.46(5), (T z ) — +0.09, and for Os 5d states (assuming 
(T z ) = 0) are fj. L = +0.016(5) and ^s = +0.096(5) giving p tot = 0.11(1). 
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ment of Os is very small and its value (not magnitude) 
grows as a function of U, passing from negative to posi- 
tive values in the case of the AMF variant of LSDA+f/. 
The magnetic dipole term on Os is predicted to be of the 
order of 10 -3 (j,b- 

The occupation of the 5/-shell of Np remains approx- 
imately 3.65, independent of the value of U and chosen 
variant of the double-counting term (FLL or AMF), al- 
though a weak decrease can be observed with increasing 
U in the FLL variant. This occupation is close to an 
expected 5/ 4 configuration (Np' ! "). consistent with the 
interpretation of the shift of the Np Mossbauer linefllj]. 
The occupation number of the 5d states of Os is also 
more-or-less insensitive to U and variant of LSDA+C/ 
potential - approximately 5.29 hd electrons are obtained 
within the Os atomic sphere (i.e., n^ s =4.71). 

Switching off the SO interaction on the Os atom does 
not seem to have a strong influence on the local charge 
distribution and magnetism of both Np and Os atoms. 
The most pronounced difference can be seen in the orbital 
magnetic moments of Np, especially in the FLL calcula- 
tions or AMF calculations with larger U . The situation 
is different when the SO interaction is switched off on Np 
atoms (while it is kept on Os atoms). In this case, we 
observe strong changes in the local electronic structures 
of both atoms. The spin moment on Np becomes sub- 
stantially larger, while the orbital moment becomes much 
more sensitive to the U value. There is also a large effect 
on the magnetic dipole term of Np, which becomes sub- 
stantially reduced. An interesting influence can be seen 
on the spin moment of the Os atoms, which become en- 
hanced by switching off the SO interaction on Np atom. 
This somewhat unusual observation can be explained by 
hybridization effects between the Np and Os electrons, 
because the spin moment grows on Np as well, providing 
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FIG. 12. (Color online) Theoretical XAS and XMCD of M 4 , 5 
edges of PuFe 2 . Calculations are presented for various values 
of the U parameter in the LSDA+C/ potential of FLL variant, 
as indicated in the legend. 



further support to the suggestion that the magnetism of 
Os is of induced nature. Switching off the SO interaction 
on both atoms almost entirely suppresses the orbital mo- 
ments and magnetic dipole terms on both atoms. In DFT 
calculations 55j the orbital and spin moments of Os as 
impurities in Fe metal are antiparallel. In NpOs2, they 
are parallel, as shown by both theory and experiment. 
This is a consequence of a complicated interplay of hy- 
bridization, Coulomb correlation, and relativistic effects, 
which arc well captured by the present calculations. 

Comparing these results to previous calculations fljj |. 
which used orbital polarization corrections (OPC) to the 
exchange-correlation potential, the individual spin and 
orbital moments on Np are of more realistic sizes in the 
present calculations. However, the OPC provides the to- 
tal moment in better agreement with experiment than 
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our LSDA+J7 calculations. It is a well-known effect that 
OPC tends to overestimate the spin and orbital poten- 
tials and, in some cases, their cancellation fortuitously 
leads to a rather good agreement with experiment, as 
found for NpOs2 [181 ] . 

For PuFc 2 , it is the FLL variant of LSDA+t/ [H| that 
provides the most satisfactory agreement with experi- 
mental spectra, see Fig. [12] The AMF LSDA+t/ cal- 
culations lead to a far too weak XMCD signal at the M5 
edge. While the FLL calculations were rather unstable 
for low values of U, at U = leV we obtained a good rep- 
resentation of the experimental XAS and XMCD spectra. 
From the theory we obtained spin moment of ^s^" Pu = 
— 1.76fj,B and orbital moment of /x^ = 2.58/is on Pu 
atoms, leading to a total moment of /i(/ t " Pu = 0.82/xs- 
The spin moment of Fe atoms is of a very similar mag- 
nitude to the Pu spin moment, namely fj,g d ~ Fc = 1.73/x_b. 
The theoretical Pu moments somewhat overestimate the 
measured ones, see Table [TtTI on the other hand the or- 
bital to spin moment ratio is in a reasonable agreement 
with experiment. 



V. DISCUSSION AND CONCLUSIONS 
A. Actinides atoms 

Earlier experimental work was done on these mate- 
rials at ANL in the 1970s |6l4l3| . and the last efforts 
in theory dates from 1990s [ljEi, SH 0, SH. In all 
cases, especially for PuFe2, where neutron experiments 
were done on a single crystal, our experimental results 
are in good agreement with previously determined quan- 
tities. XMCD has not been used previously for these 
materials and allows a more precise determination of the 
orbital spin (/is) moments and the magnetic dipole 
operator ((T z )). 

With experiments on NpOs2 we have confirmed that 
this compound has an exceptionally small total moment 
of 0.46/xs on the Np atoms. In addition, our measure- 
ments have shown that this small moment is not a result 
of a major cancellation of two large but opposing orbital 
and spin moment, but instead both of them are reduced. 
Note that DFT calculations predict for the AMF variant 
of LSDA+C/ with U = a ratio of -1.44, and for U = 1 
a ratio of -1.98, respectively. The experimental value is 
between these two, in accord with our earlier statement 
that small value of U seem to be suitable for Np-based 
Laves phases. The ratio (i 5 L f ~ Np /i/ s f ~ Np = -1.73 ±0.10, 
which is only slightly reduced from the theoretical value 
of —1.90 for IC. The reduced values of both the orbital 
and spin moments represent a particularly difficult chal- 
lenge for theory. 

The result for the branching ratios (Table [TT| gives val- 
ues of the spin-orbit operator per hole ](VF 110 } /n^-A] for 



the Np compounds of between -0.34 to -0.36, which is 
somewhat larger (by ~ 0.06) than the theoretical atomic 
value of —0.417 for IC. In good agreement with exper- 
iment, our DFT calculations suggest values between - 
0.30 and -0.38 with U between and leV. Similarly, for 
PuFe2 the value of this quantity is -0.51, that is close 
to the theoretical atomic value of -0.567 (DFT value is 
-0.50 for FLL with U=1.0eV). Since the A values are 
small, these differences could arise either from the num- 
ber of holes (10 for Np 3+ and 9 for Pu 3+ ) being too large, 
or from a change in the spin-orbit interaction caused by 
hybridization. Thus, although we are close to intermedi- 
ate coupling values throughout, there arc differences, and 
these can be ascribed to the effects of the hybridization 
with other electron states. Moreover, by comparing the 
present compounds with NP2C017 [331 ] . where the exper- 
imental value of this parameter is -0.392, such a conclu- 
sion is reinforced. The value found for NP2C017 is much 
closer to the IC theory value than found in the present 
NpX2 materials. In NP2C017 the conclusions are that the 
Np 5/ states are not strongly hybridized with the Co 3d 
states. 

Rather than discuss (T z ) we have chosen to present 
the term m m d/ '/J 5 / (recall that m m d = — 6(T 2 ), and 
pig = —2(S Z }) as this better represents the influence of 
the dipole operator, and can be more easily compared to 
theory. The term is expected to be larger for a lower sym- 
metry, as in the hexagonal symmetry found in NP2C017 
[33I ] . The value found in NpOs2 of +0.83 is slightly larger 
than the IC value (+0.56) but less than the value of +1.36 
found in NP2C017. On the basis of the values found in 
UFe2 [13] it was proposed that the expectation value of 
the magnetic dipole operator is smaller when there is 
strong hybridization and a degree of itinerancy of the ac- 
tinide 5/ electrons. There might be some evidence for 
this, especially when comparing the values in NpOs2, an 
itinerant system, with those of NP2C017, a localized sys- 
tem, but we note that the expectation value of the dipole 
operator is far from zero in NpOs2. 



B. Transition metal Os atoms 

One of the most interesting observations of the present 
study is that there is a relatively large moment induced 
on the Os atoms in NpOs2. This was not considered in 
the earlier work, and indeed demonstrates the power of 
the XMCD technique to examine sites in a compound 
other than those expected to carry magnetic moment. 
The total Os 5d moment (at an applied field of 7 T) 
is +0.09 /-is- Of this the spin moment dominates at 
+0.08 /is, and the small orbital moment is parallel to the 
spin moment, in contrast to other investigations 53, 55|. 
Thus, although hybridization is clearly present in NpOs2 
the extent and influence of both the interatomic (Np)5/ 
- (Np)5/ and the (Np)5/ - (Os)5d interactions have to 
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be understood before a full description of this material 
can be obtained. In NpOs2 the induced 5d polarization 
contributes some 30% of the total magnetism per formula 
unit, so it is by far from negligible. 

A limited number of studies have been made on the 
induced magnetization on 5d elements in ferromagnetic 
matrices. Wilhelm et al. [HI showed that for the lighter 
elements of the 5d series (specifically W) the induced spin 
and orbital moments on the ligand should be coupled an- 
tiparallcl, whereas for the heavier elements (specifically 
Ir and Pt) the coupling should be parallel. Regarding the 
induced orbital magnetic contribution, especially in the 
case of W and Ir, this contribution may be either aligned 
parallel or antiparallel to its spin induced magnetic mo- 
ment 



6lH63j. Those experimental observations of 
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a sign reversal of the orbital magnetism in transition- 
metal-based itinerant ma gne tic system were initially re- 
ported by Wilhelm et al. [53j for W in Fe/W multilayers 
and by Krishnamurthy et al. (r3^ | for Ir in Coioo-^Ir^ 
alloys. 

We compare in Fig. Ula) the XMCD spectra at the 
Os ^2,3 edges for polycrystalline Permendur doped with 
3% of Os with those for NpOs2 . In the case of this alloy, 
^d-Os = + o.50±0.005/x B /Os and the orbital contribu- 
tion is —0.045 ± 0.005/^s/Os, confirming the Tyer et al, 
[55| remark that the orbital and spin contributions should 
be anti-parallel for Os impurities in a Fe matrix. In fact, 
reference to the sum rules (Eqs. [3] & 0J shows that the 
large positive XMCD at the Os L 3 edge, is the defining 
aspect of the spectra that come from the anti-parallel na- 
ture of the spin and orbital moments in this material. In 
the case of NpOs2, however, the two moments are parallel 
in agreement with the theory presented in this paper. 

In conclusion, we have measured and calculated the 
x-ray absorption spectra and x-ray magnetic circular 
dichroism on NpAl 2 , NpFe 2 , NpOs 2 and PuFe 2 . The 
spectra calculated by density functional theory agree well 
with experimental ones for low values of Hubbard U pa- 
rameter (approximately leV). Using sum rules we have 
extracted the spin and orbital moments as well as the 
magnetic dipolc contribution of the actinide elements for 
NpOs 2 and PuFe2- We have observed a strong effect 
of surface treatment on behavior of magnetic domains, 
preventing some of the systems from reaching magnetic 
saturation even in these strong magnetic fields. Induced 
moments on Os atoms in NpOs2 were extracted and, in 
agreement with theory, the spin and orbital moments are 
aligned parallel to each other, in contrast to Os embed- 
ded in some 3d metals or alloys. 

We thank P. Voisin and S. Feite for the technical sup- 
port and the E.S.R.F. safety group, especially P. Colomp. 
We would like to thank P. Poulopoulos for providing us 
with the FeCoV:Os sample, to S. Uhle, G. Pagliosa, D. 
Bouexiere for their technical support at ITU, and to W. 
G. Stirling for encouragement in the early stages of this 



[10 

[11 
[12 

[13 
[14 

[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23; 

[24] 



K. T. Moore and G. van der Laan, Rev. Mod. Phys. 81, 
235 (2009). 

P. Santini, S. Carretta, G. Amoretti, R. Caciuffo, 
N. Magnani, and G. H. Lander, Rev. Mod. Phys. 81, 
807 (2009). 

C. Pfieiderer, Rev. Mod. Phys. 81, 1551 (2009). 

J. A. Mydosh and P. M. Oppeneer, Rev. Mod. Phys. 83, 
1301 (2011). 

H. H. Hill, in Plutonium 1970 and Other Actinides, edited 
by W. N. Miner (The Metallurgical Society of AIME, 
New York, 1970) p. 2. 

D. J. Lam and A. W. Mitchell, J. Nucl. Mat. 44, 279 
(1972). 

A. R. Harvey, Solid State Commun. 14, 1291 (1974). 
A. T. Aldred, B. D. Dunlap, D. J. Lam, and I. Nowik, 
Phys. Rev. B 10, 1011 (1974). 

A. T. Aldred, B. D. Dunlap, D. J. Lam, G. H. Lander, 
M. H. Mueller, and I. Nowik, Phys. Rev. B 11, 530 
(1975). 

A. T. Aldred, D. J. Lam, A. R. Harvey, and B. D. Dun- 
lap, Phys. Rev. B 11, 1169 (1975). 

A. T. Aldred, B. D. Dunlap, and G. H. Lander, Phys. 
Rev. B 14, 1276 (1976). 

G. H. Lander, A. T. Aldred, B. D. Dunlap, and G. K. 
Shenoy, Physica B 86-88, 152 (1977). 

A. T. Aldred, J. Magn. Magn. Mater. 10, 53 (1979). 

M. Wulff, G. H. Lander, J. Rebizant, J. C. Spirlet, 

B. Lebech, C. Broholm, and P. J. Brown, Phys. Rev. 
B 37, 5577 (1988). 

B. D. Dunlap and G. H. Lander, Phys. Rev. Lett. 33, 
1046 (1974). 

O. Eriksson, B. Johansson, M. S. S. Brooks, and H. L. 
Skriver, Phys. Rev. B 40, 9519 (1989). 
O. Eriksson, M. S. S. Brooks, and B. Johansson, Phys. 
Rev. B 41, 9087 (1990). 

O. Eriksson, B. Johansson, and M. S. S. Brooks, Phys. 
Rev. B 41, 9095 (1990). 

O. Eriksson, B. Johansson, and M. S. S. Brooks, J. Phys.: 
Condens. Matter 2, 1529 (1990). 

G. H. Lander, M. S. S. Brooks, and B. Johansson, Phys. 
Rev. B 43, 13672 (1991). 

L. Severin, M. S. S. Brooks, and B. Johansson, Phys. 
Rev. Lett. 71, 3214 (1993). 

S. P. Collins, D. Laundy, C. C. Tang, and G. van der 

Laan, J. Phys.: Condens. Matter 7, 9325 (1995). 

P. D. de Reotier, J. -P. Sanchez, A. Yaouanc, M. Fi- 

nazzi, P. Sainctavit, G. Krill, J.-P. Kappler, J. Goedkoop, 

J. Goulon, C. Goulon-Ginet, A. Rogalev, and O. Vogt, 

J. Phys.: Condens. Matter 9, 3291 (1997). 

P. D. de Reotier, J.-P. Sanchez, and A. Yaouanc, J. 



15 



Alloys Compd. 271-273, 414 (1998). 

[25] N. Kernavanois, P. D. de Reotier, A. Yaouanc, J. -P. 
Sanchez, V. Honkimaki, T. Tschentscher, J. McCarthy, 
and O. Vogt, J. Phys.: Condens. Matter 13, 9677 (2001). 

[26] M. Finazzi, P. Sainctavit, A.-M. Dias, J.-P. Kappler, 
G. Krill, J.-P. Sanchez, P. D. de Reotier, A. Yaouanc, 
A. Rogalev, and J. Goulon, Phys. Rev. B 55, 3010 
(1997). 

[27] A. Yaouanc, P. D. de Reotier, G. van der Laan, A. Hiess, 
J. Goulon, C. Neumann, P. Lejay, and N. Sato, Phys. 
Rev. B 58, 8793 (1998). 

[28] P. D. de Reotier, A. Yaouanc, G. van der Laan, N. Ker- 
navanois, J.-P. Sanchez, J. L. Smith, A. Hiess, A. Huxley, 
and A. Rogalev, Phys. Rev. B 60, 10606 (1999). 

[29] M. Kucera, J. Kunes, A.Kolmiets, M. Divis, A. V. An- 
dreev, V. Sechovsky, J. P. Kappler, and A. Rogalev, 
Phys. Rev. B 66, 144405 (2002). 

[30] F. Wilhelm, N. Jaouen, A. Rogalev, W. G. Stirling, 
R. Springell, S. W. Zochowski, A. M. Beesley, S. D. 
Brown, M. F. Thomas, G. H. Lander, S. Langridge, 
R. C. C. Ward, and M. R. Wells, Phys. Rev. B 76, 
024425 (2007). 

[31] R. Springell, F. Wilhelm, A. Rogalev, W. G. Stirling, 
R. C. C. Ward, M. R. Wells, S. Langridge, S. W. Zo- 
chowski, and G. H. Lander, Phys. Rev. B 77, 064423 
(2008). 

[32] T. Okane, T. Ohkochi, T. Inami, Y. Takeda, S. i. Fu- 
jimori, N. Kawamura, M. Suzuki, S. Tsutsui, H. Yam- 
agami, A. Fujimori, A. Tanaka, D. Aoki, Y. Homma, 
Y. Shiokawa, E. Yamamoto, Y. Haga, A. Nakamura, and 
Y. Onuki, Phys. Rev. B 80, 104419 (2009). 

[33] I. Halevy, A. Hen, I. Orion, E. Colineau, R. Eloirdi, J.- 
C. Griveau, P. Gaczyski, F. Wilhelm, A. Rogalev, J.-P. 
Sanchez, M. L. Winterrose, N. Magnani, A. B. Shick, 
and R. Caciuffo, Phys. Rev. B 85, 014434 (2012). 

[34] T. Shishidou, T. Oguchi, and T. Jo, Phys. Rev. B 59, 
6813 (1999). 

[35] V. N. Antonov, B. N. Harmon, and A. N. Yaresko, Phys. 

Rev. B 68, 214424 (2003). 
[36] V. N. Antonov, B. N. Harmon, O. V. Andryushchenko, 

L. V. Bekenev, and A. N. Yaresko, Phys. Rev. B 68, 

214425 (2003). 

[37] A. N. Yaresko, V. N. Antonov, and B. N. Harmon, Phys. 

Rev. B 68, 214426 (2003). 
[38] V. N. Antonov, B. N. Harmon, O. V. Andryushchenko, 

L. V. Bekenev, and A. N. Yaresko, Low Temp. Phys. 30, 

305 (2004). 

[39] J. Kunes, P. Novak, M. Divis, and P. M. Oppeneer, Phys. 

Rev. B 63, 205111 (2001). 
[40] B. T. Thole, P. Carra, F. Sette, and G. van der Laan, 

Phys. Rev. Lett. 68, 1943 (1992). 
[41] P. Carra, Synchrotron Radiation News 5, 21 (1992). 
[42] P. Carra, B. T. Thole, M. Altarelli, and X. Wang, Phys. 

Rev. Lett. 70, 694 (1993). 
[43] G. van der Laan and B. T. Thole, Phys. Rev. B 53, 14458 

(1996). 



[44] T. Okane, Y. Takeda, J. Okamoto, K. Mamiya, 
T. Ohkochi, S. i. Fujimori, T. Saitoh, H. Yamagami, 
A. Fujimori, A. Ochiai, and A. Tanaka, J. Phys. Soc. 
Jpn 77, 024706 (2008). 

[45] A. Rogalev, J. Goulon, C. Goulon- Ginet, and C. Mal- 
grange, Lect. Notes Phys. 565, 61 (2001). 

[46] J. Goulon, C. Goulon-Ginet, R. Cortes, and J. M. 
Dubois, J. Phys. 43, 539 (1982). 

[47] L. Troger, D. Arvanitis, K. Baberschke, H. Michaelis, 
U. Grimm, and E. Zschech, Phys. Rev. B 46, 3283 

(1992) . 

[48] P. Pfalzer, J.-P. Urbach, M. Klemm, S. Horn, M. L. den- 
Boer, A. I. Frenkel, and J. P. Kirkland, Phys. Rev. B 
60, 9335 (1999). 

[49] M. J. Berger, J. H. Hubbell, S. M. Seltzer, J. Chang, 
J. S. Coursey, R. Sukumar, D. S. Zucker, and K. Olsen, 
(2010), [http://www.nist.gov/pml/data/xcom]. 

[50] The edge jump ratio of the M5 to M4 edges can be de- 
duced from the entries in these tables for all actinides. 
They are all 1.57±0.05, except for Pu, which is listed as 
1.87. We believe this value is incorrect and we have used 
1.57 for our data analysis. 

[51] G. van der Laan, K. T. Moore, J. G. Tobin, B. W. Chung, 
M. A. Wall, and A. J. Schwartz, Phys. Rev. Lett. 93, 
097401 (2004). 

[52] A. J. Freeman, J. P. Desclaux, G. H. Lander, and 
J. Faber, Phys. Rev. B 13, 1168 (1976). 

[53] F. Wilhelm, P. Poulopoulos, H. Wende, A. Scherz, 
K. Baberschke, M. Angelakeris, N. K. Flevaris, and 
A. Rogalev, Phys. Rev. Lett. 87, 207202 (2001). 

[54] This was deduced from the changes in white lines inten- 
sities in NpOs2 with respect to Os impurities in the per- 
mendur reference sample, in which the estimated number 
of 5d holes is n| d =3.70. 

[55] R. Tyer, G. van der Laan, W. M. Temmerman, Z. Szotek, 
and H. Ebert, Phys. Rev. Lett. 90, 129701 (2003). 

[56] P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, 
and J. Luitz, WIEN2k, An Augmented Plane Wave + Lo- 
cal Orbitals Program for Calculating Crystal Properties, 
Techn. Universitat Wien, Austria (2001). 

[57] M. T. Czyzyk and G. A. Sawatzky, Phys. Rev. B 49, 
14211 (1994). 

[58] V. I. Anisimov, I. V. Solovyev, M. A. Korotin, M. T. 
Czyzyk, and G. A. Sawatzky, Phys. Rev. B 48, 16929 

(1993) . 

[59] J. Kunes, P. Novak, R. Schmid, P. Blaha, and 
K. Schwarz, Phys. Rev. B 64, 153102 (2001). 

[60] J. Kunes, P. Novak, R. Schmid, P. Blaha, and 
K. Schwarz, Phys. Rev. B 64, 153102 (2001). 

[61] G. Schiitz, M. Kniille, and H. Ebert, Phys. Scr. T49A, 
302 (1993). 

[62] V. V. Krishnamurthy, D. J. Singh, N. Kawamura, 
M. Suzuki, and T. Ishikawa, Phys. Rev. B 74, 064411 
(2006). 

[63] M. A. Laguna-Marco, D. Haskel, N. Souza-Neto, J. C. 
Lang, V. V. Krishnamurthy, S. Chikara, G. Cao, and 
M. van Veenendaal, Phys. Rev. Lett. 105, 216407 (2010). 



